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Nicolis and Piazza have recently pointed out the existence of Nambu-Goldstone-like excitations in relativistic 
systems at finite density, whose gap is exactly determined by the chemical potential and the symmetry algebra. 
We show that the phenomenon is much more general than anticipated and demonstrate the presence of such 
modes in a number of systems from (anti)ferromagnets in magnetic field to superfluid phases of quantum chro- 
modynamics. Furthermore, we prove a counting rule for these massive Nambu-Goldstone bosons and construct 
a low-energy effective Lagrangian that captures their dynamics. 
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Introduction. — Trying to understand collective behavior of 
matter in nonlinear many-body systems is a challenge com- 
mon to many areas of physics. At long distances and low 
temperatures, excitations with vanishing or small gap (mass) 
dominate the dynamics. The concept of spontaneous sym- 
metry breaking has been cruciai for its understanding, as 
it unambiguously predicts existence of gapless excitations — 
the Nambu-Goldstone bosons (NGBs) — such as phonons or 
magnons. For nearly five decades, however, their correct 
counting and dispersion relations eluded consistent under- 
standing. Recently, the present authors developed a unified 
framework to determine the number and dispersion relations 
of NGBs |HI2l, including their redundancies fl. 

Cases where exact statements can be made about gapped 
modes are rare though. Kohn's theorem states that a gas of 
charged particles with Galilean invariance, when exposed to a 
uniform magnetic field, sustains a collective mode with the cy- 
clotron gap [4]. Moreover, some soliton solutions to nonlinear 
equations saturate Bogomol'nyi-Prasad-Sommerfieldbounds, 
allowing their energies to be determined based on symmetry 
alone [5], albeit with limited applicability to observable sys- 
tems. NGBs perturbed by explicit symmetry breaking effects 
acquire small gaps and are usually called pseudo-NGBs [6]. 
Yet their gaps in general can be computed only approximately. 

Recently, Nicolis and Piazza J3l pointed out that the gaps 
of pseudo-NGBs can be determined in special circumstances. 
Considering Lorentz-invariant systems perturbed only by a 
chemical potential whose charge operator is spontaneously 
broken, they showed that the mass of the associated pseudo- 
NGB can be computed exactly and is free of radiative correc- 
tions. We will cali such state massive NGB (mNGB). In the 
present Letter, we show that mNGBs appear in a much broader 
class of systems; the theory need not be Lorentz-invariant, or 
the chemical potential operator spontaneously broken. We 
prò vide a counting rule for the number of mNGBs and con- 
struct an effective Lagrangian description for them. 



General argument. — Consider a many-body system speci- 
fied by the Hamiltonian H with an internai symmetry group 
G. In order to describe states with finite charge density, it 
is customary to introduce a chemical potential fi by H = 
H — fiQ, where Q is one of the generators of G. The vacuum 
|0) is defined as the eigenstate of H with the lowest eigen- 
value. Without lack of generality, we can take H\0) = 0. 
Since the generators Qi of the Lie group G commute with the 
Hamiltonian H, they are ali time-independent in the Heisen- 
berg picture defined by H fl, 
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where jf(x) are the corresponding locai charge densities. 

When spontaneously broken, generators of the symmetry 
group G of the full Hamiltonian H give rise to standard mass- 
less NGBs. On the other hand, spontaneously broken genera- 
tors that do not commute with Q create pseudo-NGBs whose 
gaps are proportional to fi. The novelty of the observation 
made by Nicolis and Piazza |7] is that the proportionality con- 
stants can be computed exactly by group theory. 

By the standard Cartan decomposition, explicitly broken 
generators can be split into pairs Q± a — the roots — such that 



[Q,Q±o] = ±q<jQ±o, 



(2) 



where Q± a are some complex linear combinations of explic- 
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Let us now focus 



on the quantity X a = (0\[Q +<J (t), j" CT (0)]|0), which is man- 
ifestly time-independent. Using Eq. (Q~|), inserting a complete 
set of eigenstates \n,p) of momentum P and energy H, and 
carrying out integration over space, we obtain 
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Provided that pq a > 0, time-independence of the left-hand 
side and E n (0) > require (0|j°_ (T (0)|n, 0) = for each n. 
If A(j is zero, implying (0|j'9 CT (0)|rì, 0) = for each n as well, 
Q aR = (Q +a + Q- a )/V2 and Q aI = (Q +<7 - Q_ CT )/\/2i 
cannot be spontaneously broken. Namely, there is no locai 
fìeld $(x) suchthat (0\[Q„ RJ (t), $(0)]|0) ^ 0. 

On the other hand, if X a ^ 0, Q a R and Q ff / are broken 
spontaneously and there must be a state \n, 0) with mass 



H\n,0) = £ n (0)|n,0) = W(r |n,0} 



(4) 



suchthat (0|j%(0)|n,0) ^ and (0|j° CT (0)|n, 0) = 0. This 
is the mNGB associated with the pair Q± a . 

Our derivation clarifìes several points on mNGBs. First, it 
is clear that the assumption of the underlying dynamics being 
Lorentz-invariant U7|] can be dropped. Also, A CT always plays 
the role of the order parameter for charges Q a R.i- Finally, 
note that (0\j°_ CT (0)|n, 0) = for ali n means Q +cr |0> = 0, 
while if A,j is nonzero, Q-er|0) 7^ 0. This observation leads 
to a simpler, albeit less rigorous, understanding of mNGBs. 
Eq. d2J gives [H, Q± a ] = Tpq<jQ±a which implies that the 
state Q-<j|0) has energy pq a . As there cannot be a state with 
energy lower than the vacuum, Q +tT \0) has to vanish. Our 
argument is reminiscent of the proof of Kohn's theorem [4], 
allowing for a unifìed comprehension of the two phenomena. 

Number of mNGBs. — For a proper understanding of the 
low-energy dynamics of the system, it is important to know 
the number and dispersion relations of NGBs. Denotine the 
broken generators of G as Q a , the former is given by I1II2J], 

«-NGB = n A + n B , n A — n BC - rank p, n B = ìrankp, 

(5) 

where n B c is the number of broken generators and 
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[Qa,Qb]\0), 



(6) 



il being the spatial volume. The type-A and B NGBs gener- 
ically have linear and quadratic dispersion relations and cor- 
respond to type-I and II in the Nielsen-Chadha theorem J9|, 
even though this is not always the case IlOll . Each type-B 
NGB is described by a canonically conjugate pair of broken 
generators Q a and Qi, with nonzero p a {, and hence two broken 
symmetries count as one degree of freedom, whereas type-A 
NGBs are stand-alone like in the originai Goldstone theorem. 

Here we address the question of counting the mNGBs 11 ili . 
Namely, we show that their number is given by 



rimNGB = -(rankp - rankp), 



(7) 



where the matrix p is defined analogously to Eq. © for ali 
generators of G instead of just G. To that end, we have 
to further specify the structure of the Lie algebra. First, 
let us choose the maximal number of mutually commut- 
ing generators of G, including Q itself, to form the Car- 
tan subalgebra. By a proper choice of the vacuum |0), 
we can always ensure that these are the only generators 



of G that can possibly have a nonzero vacuum expecta- 
tion value ti]. This alone does not preclude the possibil- 
ity that explicitly broken generators acquire expectation val- 
ues. However, since ±pq a (0\Q± a \0) = (0\\pQ,Q± a ]\0) = 
(0\[Q± a ,H]\0) = thanksto H\0) = 0, (0|Q ±(T |0) must 
vanish for any nonzero q a . If we arrange the generators as 

Qi = (QiR,Qu,---,QmR,Qmi,Qi,---,Q d i m G)> where 
ni = (dimG — dimG)/2, the matrix p becomes block- 
diagonal, p = diag(iT2Ai, . . . , ÌT2X m , p), t~2 being the second 
Pauli matrix. Thus, \ (rank p — rank p) counts the number of 
nonzero Acr's. Assuming that there is at most one mNGB for 
each pair of Q± a , this proves our counting mie d7). In the fol- 
lowing, we provide examples of mNGBs, demonstrating the 
validi ty of Eq. (Q in physically interesting systems 111 211 . 

Ferromagnet. — The Hamiltonian of a ferromagnet enjoys 
the internai G = 0(3) symmetry group of spin rotations. In 
the ground state, individuai spins are aligned, breaking this 
symmetry down to its 0(2) subgroup. The two broken gen- 
erators give rise to a single type-B NGB with a quadratic dis- 
persion relation at low momentum — the magnon J2, IH . 

Consider now switching on a uniform magnetic fìeld B ori- 
ented in the z-direction. This amounts to breaking the sym- 
metry explicitly to G = 0(2) by adding to the Hamiltonian 
—p m BS z (p m B > 0), where S is the total spin operator and 
p m is the magnetic moment. This term can be viewed as a 
chemical potential p = p m B for the generator Q = S z . 
Given that [S Z ,S±] = ±S± where S± = {S x ± iS v )/y/2, 
S- must excite a mNGB of gap p, which is just the magnon 
with energy lifted by the magnetic fìeld 11 1 311 . The operator S+ 
annihilates the ground state. Both these assertions are easy to 
understand from the fact that the vacuum corresponds to the 
state with maximum spin in the direction of the magnetic fìeld, 
and the magnon to an excitation caused by flipping one of the 
spins. Note that the counting mie (Q predicts the correct num- 
ber of mNGBs, that is, T^mNGB = = 1. 

Antiferromagnet. — In the absence of a magnetic fìeld, as- 
sume the spins are oriented alternately along the z-axis; G — 
0(3) is broken to 0(2) just like in a ferromagnet. In this case 
there are two type-A NGBs, one for each broken generator. 

Applying a magnetic fìeld along the z-axis leads to an insta- 
bility as the NGBs attempt to acquire masses ±p = ±p m B. 
The ground state rearranges with altemating spins pointing 
in an orthogonal direction instead, say along the x-axis. Then 
Q = S z is a spontaneously broken generator which commutes 
with H and creates a gapless type-A NGB. On the other hand, 
the other spontaneously broken generator S y is explicitly bro- 
ken at the same time, creating a mNGB with gap p lfl4ll . 
The magnetic fìeld induces a small magnetization along the 
z-axis, and hence p xy — (Q\[S X , S y ]\Q) ^ 0. Consequently, 
«mNGB = = 1, consistent with Eq. (|7]). 

Relativistic Bose-Einstein condensation. — As an explicit 
example where p 7^ 0, consider a theory of a complex scalar 
doublet 4> with a global G = SU(2) x U(l) symmetry, 

Jgf = D^D^cj)- M 2 tf<t>- A(0 f 0) 2 , (8) 

where Dqc/) = (do — ip)(f>. This model features a relativis- 
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tic Bose-Einstein condensation (BEC) phase for fi > M, in 
which the symmetry is spontaneously broken to a U(l)' sub- 
group. The three broken generatore produce two NGBs, one 
type-A and one type-B 11511 . consistent with Eq. ((5]), since one 
of the SU(2) charges develops nonzero density in the ground 
state, hence rank p = 2. 

The type-B NGB in this model has an "antiparticle", car- 
rying opposite charge. Its mass equals 2/i and does not re- 
ceive radiative corrections 11 1611 . To see why, note that when 
fi = 0, the Lagrangian enjoys an extended internai symmetry, 
G = SO(4) ~ SU(2) L x SU(2) R . This is most easily seen 
by defining a 2 x 2 matrix $ = (</>, ÌT20*), which transforms 
under G as $ — > Ul$U^. Denote the generatore of G as L and 
R, respectively; they are both given by a half of the Pauli ma- 
trices. The SU(2) rotations of the doublet now correspond 
to SU(2)l; the U(l) phase transformations are generated by 
2i?3. The remaining two generatore of SU(2)r are explicitly 
broken by the chemical potential fi. In the BEC phase, the 
condensate can be chosen as (0|$|0) ~ 1 so that they are also 
broken spontaneously. Since the R± generators of SU(2)r 
satisfy the commutation relation [2i?3, R±] = ±2R±, Eq. (01 
implies that i?_ creates a mNGB with mass 2 fi, in agreement 
with the explicit calculation. Indeed, mmUGB = -^p = 1- 
This example obviously admits a generalization to a large 
class of relativistic linear sigma models with chemical po- 
tential 11171 [l8Q , the key ingredient being an extended global 
symmetry when the chemical potential is set to zero. 

QCD-like theories. — Quantum ChromoDynamics (QCD) 
with two degenerate quark flavors possesses an approximate 
global SU(2)l x SU(2)r chiral symmetry. A nonzero quark 
mass breaks this explicitly to the G = SU(2)v subgroup gen- 
erated by V = R + L. The chiral condensate in the QCD 
vacuum breaks the symmetry spontaneously in the same way, 
resulting in three pseudo-NGBs in the spectrum: the pions. 

Nonzero chemical potential, fii, for V3 breaks the exact 
symmetry G further to its G — U(l)i subgroup, generated 
by V3. While the mass of the neutral pion is insensitive to 
(il, the masses of the charged pions become m n ± fi\. Conse- 
quently, once fii > m„ the positively charged pion undergoes 
BEC, breaking the residuai G symmetry spontaneously 11911 . 
Therefore, the spectrum in the pion BEC phase exhibits one 
true, type-A NGB. However, the ground state has a nonzero 
isospin density, (0|V 3 |0) = -i{Q\\Vi, V 2 }\0), and Eq. @ im- 
plies that there is also one mNGB. Such a state has indeed 
been found using effective field theory [20] as well as vari- 
ous model approaches 11 8l felli and can be identifìed with the 
neutral pion in the superfluid medium. As opposed to these 
approximate calculations, the result of Ref. [7] nevertheless 
ensures that its mass is exactly equal to fi\. This follows from 
the commutation relation [V3, V±] = ±V±. 

In the limit of massless quarks, the full symmetry be- 
comes G = SU(2)l x SU(2)r; isospin chemical potential 
breaks this explicitly to G = U(1)l X U(1)r. Pion conden- 
sate now develops at any nonzero chemical potential, break- 
ing G spontaneously to U(l). Thus, there is one type-A 



NG boson in the spectrum. Moreover, given the commuta- 
tore [V3,i?±] = ±R± and [V3,£±] = ±L±, we find that 
(0|F 3 |0) = -2»<0|[iJi,ik]|0) = -2i(0|[ii,i a ]|0) ^ 0, as 
a result of which there are = 2 mNGBs according to 
Eq. (0. This is consistent with explicit calculations; the addi- 
tional mNGB has the quantum numbers of the a meson. 

A related theory that has attracted attention in connection 
to lattice Monte-Carlo simulations at nonzero density is two- 
color QCD. Thanks to extended global symmetry, its two- 
flavor version has altogether fìve pseudo-NGBs in the vac- 
uum: three pions and an (anti)diquark 02211 . For baryon chem- 
ical potentials /ìb exceeding m^, the diquark will thus con- 
dense. Using the same argument as above, one can show that 
in this BEC phase, the pions form a triplet of mNGBs with 
mass exactly equal to fi^. Such modes have been observed 
in analytic calculations 122, 23tl as well as on the lattice Ì24ll . 
In the limit of massless quarks, the enlarged symmetry again 
leads to the presence of an extra mNGB with the quantum 
numbers of the a meson. Similar conclusions can be reached 



for an arbitrary even number of flavors [22 1. 

Effective Lagrangian formalism. — For energies, including 
the chemical potential fi, well below the characteristic scale 
of spontaneous symmetry breaking, the dynamics of both 
NGBs and mNGBs can be described using effective field the- 
ory 1 1 3l 12511 . Assuming spatial rotational invariance and using 
the fact that the chemical potential can be regarded as a Con- 
stant temporal gauge field [26], Aq — fi 1 (i — 1, . . . , dim G), 
the leading-order effective Lagrangian reads fi 1 311 



^eff =c a (7r)7r a + e^ir)^ + l -g ah {n)D t w a D t TT b 
-\g ab (K)VTT a -VTT b . 



(9) 



The NG fields n a (a = 1, . . . , dim G/H) live on the coset 
G/H of the broken symmetry, H being the unbroken sub- 
group at zero chemical potential; g a b(^) and g a b(^) are both 
G-invariant metrics on the coset. Under an infinitesimal sym- 
metry transformation defined by a set of parameters e 1 , the 
coset fields change as Sn a — e l hf(Tr), where hf(ir) are the 
Killing vectors of the metrics. The covariant derivative is 
D t ir a = ii a - fi l h?(ir). 

Explicit expressions can be obtained using the formalism 
developed in Ref. lfe~7ll . Denoting now the broken group gen- 
erators as T a and the unbroken generators as T p , we represent 
the coset element by U(tt) = é lTa7T and define the Maurer- 
Cartan form as u; a (7r) = Titu^ir) = —iU(ir)ì-£^;U(Tr). As- 
suming that G is compact and that the generators are normal- 
ized according to Tr(XÌTj) = XSij, we obtain 112811 . 



?a&M = 5c<i(0K(7rK(7r), e 4 (7r) = i/f (7^(0), 



(10) 



where vì(%) = iTr[C/(7r) t T i C/(7r)r :; ]. For consistency with 
the G-invariance of the action, the effective couplings ej(0) 
and g ab (0) must satisfy f./e^O) = and f pa c g cb {0) + 
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/ pb c g ac (0) = 0, where f i - k are the structure constants of the 
symmetry group. Similar expressions hold for g. 

Let us emphasize that while the effective Lagrangian (0 
was derived by Leutwyler nearly two decades ago [131, he 
only obtained a set of differential equations determining the 
effective couplings c a (n), e^w), g a b{^), 9ab(^) implicitly. To 
the best of our knowledge, this is the first time that a general 
solution has been found, maintaining the full nonlinear depen- 
dence of the couplings on the coset fields 7r a . Using Eq. ( fTot . 
the effective Lagrangian is now completely fìxed by the values 
of ffab(O) and g a i,(0), encoding decay constants of the NGBs, 
and of ei(0), expressing charge densities in the ground state. 

The effect of the chemical potential on the ground state and 
spectrum can be completely deduced from these expressions. 
First of ali, the Lagrangian © acquires a potential 

V(n) - -^(7r) e< (0) - \gab(0)v a (n)v b (Tr), (11) 

where v % (tt) = fjPi/Un). As observed above, turning on the 
chemical potential can lead to an instability of the vacuum. 
The true ground state is then found by minimizing the po- 
tential. In particular, when ei(0) = for ali i, chemical 
potentials for unbroken generators T p generate the potential 
-Ì5cd(0)(/i p / a /^ a )(M ,T / 6CT <i ^) + O(T 4 ),whichis negative 
definite so that NGBs condense. This is what happens when 
an antiferromagnet is subjected to a magnetic fìeld along the 
z-axis: BEC of NGBs leads to a (canted) Néel order in the xy 
piane. On the other hand, when (0) ^ for some i, BEC 
of NGBs always results in a new ground state where /i 1 is par- 
allel to ei(0), at least for an infìnitesimal //. For example, in 
the case of ferromagnet where e z (Q) = (0\S z \0)/il > dom- 
inates, a magnetic fìeld in ir-direction causes rearrangement 
of the magnetization into .x-direction, while a magnetic fìeld 
in z-direction further stabilizes the ground state. 

The case of ferromagnet admits a straightforward general- 
ization. Consider a system which features: (i) an unbroken 
U(l) symmetry; (ii) a type-B NGB charged under it; (iii) no 
type-A NGBs charged under it. Then, introducing a chemi- 
cal potential for the unbroken subgroup with a suitably cho- 
sen sign will lift the energy of the type-B NGB, making it a 
mNGB. No instability occurs since a type-B NGB is described 
by a pair of canonically conjugate NG fields; there is no gap- 
less mode charged oppositely under the unbroken U(l). 

Provided that (0) = for ali i, simple results can be ob- 
tained also in the case that chemical potentials are introduced 
for the broken generators T a . Assuming that the coset G/H 
is a symmetric space and NGBs transforms irreducibly under 
the unbroken subgroup H , the first-order derivative term starts 
from 0{ir 3 ). The potential term reads \\i c >° ( j ' C J '] ' db ' W \ 
explaining the mass of the mNGBs. Of course, the applicabil- 
ity of our effective Lagrangian is not restricted to this partic- 
ular case; we checked that ali examples we mentioned above 
can be explained by using our formalism. 

Quasìparticle nature of mNGBs. — What happens to the 
mNGB when fi exceeds the symmetry-breaking scale? The 
general argument asserts just the existence of a state which, in 



the long-wavelength limit, is an exact eigenstate of the Hamil- 
tonian. True NGBs are actually well-defined quasiparticles 
even at finite wavelength, as the ratio of their width and their 
energy goes to zero in the zero-momentum limit [16]. 

First, note that in ali the examples discussed above, the 
mNGB is stable at low enough chemical potential. The ex- 
amples where Q is not spontaneously broken are trivial as the 
mNGB is then the lightest mode carrying its charge. Even in 
other cases, the mNGB is protected against decay by an un- 
broken symmetry. Indeed, the mNGB in the SU(2) x U(l) 
linear sigma model (0 is stable due to conservation of the un- 
broken U(l)' charge. The pion superfluid phase of QCD with 
isospin chemical potential has an unbroken discrete symme- 
try: a combination of parity with a U(l)i rotation by the angle 
ir; the neutral pion is the only NGB mode odd under this trans- 
formation. Finally, in the diquark BEC phase of two-color 
QCD with baryon chemical potential, the pion triplet is stable 
as a consequence of isospin and parity conservation. 

However, the mNGB is not expected to remain stable for 
arbitrarily high chemical potential. Consider the case of QCD 
with isospin chemical potential as an example. At very high 
/ìi, the system behaves as a weakly coupled superfluid. The 
quarks are gapped, but although their gap A increases with 
/il, the ratio A//ìi tends to zero with increasing /ìi 11911 . The 
mNGB will thus dissolve into a pair of quarks. This is in no 
contradiction with the general result, only its implications are 
trivial at high chemical potential: the mNGB excitation is real- 
ized as a state in a two- or many-body continuum. In general, 
the mNGB mode is expected to be a well-defined quasìparti- 
cle only for chemical potentials well below the characteristic 
scale of spontaneous symmetry breaking. 

Conclusione. — We have generalized the argument of 
Ref. J3l and shown that a large class of quantum many-body 
systems features mNGBs whose mass is exactly determined 
by the chemical potential and the symmetry algebra. Since the 
chemical potential acts as a Constant background gauge fìeld, 
this can be understood as a consequence of gauge invariance, 
encoded in the effective Lagrangian ((9). Moreover, we have 
derived a general expression for the number of mNGBs. 

Of course, the calculation of the mass and the discussion of 
stability presented here is merely the first step towards a more 
detailed analysis of the properties of mNGBs. It would be cer- 
tainly interesting to investigate which of the other hallmarks 
of NGBs are preserved in presence of the chemical potential 
explicitly breaking part of the symmetry. These include for in- 
stance the saturation of Ward-Takahashi identities for correla- 
tion functions of the broken currents and low-energy theorems 
for the scattering ampli tudes involving one or more NG parti - 
cles. We expect that in ali these future efforts, the underlying 
gauge invariance will play a centrai role. 
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